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Abstract: The Crab Pulsar Wind Nebula is the best studied source of 7-ray astrophysics. The contribution of the various 
soft radiation fields to the Inverse Compton component of its high energy emission, the strenght of the internal magnetic 
field and the maximum energies reached by primary electrons are however still matter of study. 

The MAGIC stereoscopic system recorded almost 50 hours of Crab Nebula data in the last two years, between October 
2009 and April 201 1. Analysis of this data sample using the latest improvements in the MAGIC stereo software provided 
an unprecedented differential energy spectrum spanning three decades in energy, from 50 GeV up to 45 TeV. At low ener- 
gies, the MAGIC results, combined with the Fermi/LAT data, yield a precise measurement of the Inverse Compton peak. 
In addition, we present light curves of the Crab Nebula at different time scales, including a measurement simultaneous 
to one of the Crab Nebula flares recently detected by both Fermi/LAT and AGILE. Using the MAGIC spectrum together 
with multiwavelength data, we discuss the implications for the modeling of the Crab Nebula. 
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1 Introduction 

The pulsar wind nebula (PWN) associated with the Crab 
Pulsar is a leftover of a supernova explosion recorded in 
1054 [Ij. Located at the center of the nebula, the pul- 
sar continuously releases its energy primarily in the form 
of a highly relativistic and magnetized wind of particles, 
mainly electrons and positrons. This pulsar wind termi- 
nates in a standing shock where particles are accelerated up 
to ultra-relativistic energies, and their pitch angles are ran- 
domized. The outflow interacts with the surrounding mag- 
netic and photon fields creating the PWN. The nebula emits 
synchrotron radiation which is observed from radio fre- 
quencies up to soft 7-rays. This emission is well-described 
in terms of the magnetohydrodynamic model (MHD) 1 2 1 . 
At higher energies (above 500 MeV), the overall emission 
is dominated by the Inverse Compton (IC) scattering by 
the electrons from the pulsar, predominantly on the syn- 
chrotron photons 1 3 1 . 

The Crab Nebula is one of the best-studied non-stellar ob- 
jects in the sky. Due to its bright glow at almost all wave- 
lengths, precise measurements can be provided by many 
different kinds of instruments, allowing for a detailed ex- 
amination of its physics. 

The IC emission from the Crab Nebula was detected for 



the first time above 700 GeV by the pioneering Whipple 
telescope in 1989 |4|. Since then, the imaging Cherenkov 
technique was successfully used to extend the Crab Neb- 
ula differential energy spectrum from a hundred GeV up to 
a hundred TeV. Nevertheless, the low-energy part of this 
component, below 100 GeV, where the IC peak resides, 
could be observed only in the last years. At low energies, 
the Fermi satellite filled the gap between few and hundred 
GeV |5|. In the meanwhile, the construction of imaging 
atmospheric Cherenkov telescopes (lACTs) with larger re- 
flective surface allowed to lower the energy threshold of the 
ground-based telescopes below 100 GeV. The stand-alone 
first MAGIC telescope had already shown a hardening of 
the spectrum below a few hundreds GeV |6 |. With an en- 
ergy threshold of 80 GeV, MAGIC provided the only exper- 
imental data points overlapping with the Fermi/LAI ones 
up to now. The goal of this work is to measure the Crab 
Nebula differential energy spectrum with a higher statisti- 
cal precision and down to lowest possible energy threshold 
by using the MAGIC stereoscopic system. 

In addition, because of its apparent overall flux steadiness, 
the Crab Nebula has been considered as standard candle at 
almost all wavelengths. It has been used to cross-calibrate 
X-ray and 7-ray telescopes, to check the instrument perfor- 
mance over time, and to provide units for the emission of 
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other astrophysical objects. However, this notion surpris- 
ingly changed in September 2010, when both AGILE and 
FermilhPil satellites observed an enhancement of the 7-ray 
flux above 100 MeV by a factor 2 |7, 8|. Spectral analysis 
revealed that the flare had a synchrotron origin, and a spec- 
trum which is harder than the average one. Optical and 
X-ray images taken right after the flare showed a few neb- 
ular, brightened features. One of them is the "anvil", the 
knot which lies in the inner nebula in the projected inward 
extension of the jet, and it is considered a primary site for 
particle acceleration. If the 7-ray flare luminosity suggests 
that the production region is close to the pulsar, the short 
flare rising time favors a compact emission region of size 
smaller than 10^^ cm. Thus, the "anvil" feature could be 
an excellent flare site candidate I?). A similar episode was 
observed in April 2011 by the same teams (ATel ^3276). 
It lasted for two weeks, with the peak occurring on April 
14, 2011. 

The MAGIC telescopes recorded some data simultaneously 
to both high-energy 7-ray flares in September 2010 and 
April 2011. The results of these observations will also be 
presented in this paper. 

2 The MAGIC telescopes 

MAGIC consists of two 17 m diameter lACTs located in 
the Canary island of La Palma, Spain (2200 m above sea 
level). It became a stereoscopic system in Autumn 2009. 
The stereoscopic observation mode led to a significant im- 
provement in the performance of the instrument |9||10|. 
The current sensitivity of the array yields 5a significance 
detections above 250 GeV of fluxes as low as 0.8% of the 
Crab Nebula flux in 50 hr. The energy threshold of the 
analysis lowered down to 50 GeV for low zenith angle ob- 
servations. 

For this work we considered stereoscopic observations of 
the Crab Nebula carried out during the last two years, be- 
tween October 2009 and March 2011. These observations 
are performed in wobble mode, at zenith angles between 5° 
and 50°. Data affected by hardware problems, bad atmo- 
spheric conditions, or displaying unusual background rates 
were rejected in order to ensure a stable performance. 48.7 
hours of effective time are left after data quality selection. 

The analysis was performed by using the tools of the 
MAGIC analysis software ifTDl . Each telescope records 
only the events selected by the hardware stereo trigger. 
The so-called image cleaning procedure selects the pixels 
which have significant signal. The obtained reconstructed 
image is then parametrized with a few simple quantities. In 
particular, we use a new algorithm, the so-called sum im- 
age cleaning to improve the analysis sensitivity at energies 
below 100 GeV. Similarly to the standard algorithm it uses 
the concept of core and boundary pixels, defined according 
to their charge levels. The addition of new constraints on 
the charge level of compact group of neighboring core pix- 
els, within a very short time spread, further allows to lower 
the charge thresholds. This is an important aspect in the 



reconstruction of low-energy shower images which contain 
only few pixels. The sum image cleaning performs as well 
as the standard one for energies above 150 GeV, whereas 
below such energies it improves the sensitivity by a 15%. 

Afterwards, pairs of images from the same stereo events 
are combined and the shower direction is determined as 
the average of the corresponding single-telescope direc- 
tions LIOJ . The background rejection relies on the defi- 
nition of the multi-variable variable hadronness, which is 
computed by means of a Random Forest (RF) algorithm. 
RF uses as input a small set of image parameters from both 
telescopes, together with the information about the shower 
location provided by stereoscopy. The 7-ray signal is esti- 
mated through the distribution of the squared angular dis- 
tance (0^) between the reconstructed and the catalog source 
position. The energy of each event is estimated by using 
look-up tables created from Monte Carlo simulated 7-ray 
events. For the computation of the differential energy spec- 
trum, the 7-ray signals in each energy bin are determined 
by selecting a soft hadronness cut which retains 90% of the 
7-ray events. 

The data sample was divided in four zenith angle ranges to 
account for corresponding variations in the image param- 
eters. The matrices for the background rejection obtained 
through the RF were computed for each sub-sample sepa- 
rately. The four independent analyses were combined later 
at the last stage of the analysis. In Figure [T] we show the 
collection area obtained for the four data sub-samples. 
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Figure 1 : Collection area obtained for the data sub-samples 
in the four different ranges of zenith angle. 



3 The differential energy spectrum 

Figure |2] shows an unprecedented differential energy spec- 
trum of the Crab Nebula which spans three orders of mag- 
nitude in energy, from 50 GeV up to 45 TeV In order to 
correct for the finite energy resolution and energy estima- 
tion bias it was unfolded with Bertero's method 1 13 1. Three 
other unfolding methods were considered, and all of them 
gave compatible results within the statistical errors. The 
yellow-shadowed area indicates the systematic uncertain- 
ties on the flux normalization and the spectral slope. The 
first one is estimated to be 15% of the flux value, whereas 
the second one produces an absolute uncertainty on the 
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Figure 2: Differential energy spectrum of the Crab Nebula 
obtained with data recorded by the MAGIC stereoscopic 
system. 
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Figure 3: Spectral energy distribution of the Crab Nebula 
obtained with the MAGIC telescopes, together with the re- 
sults from previous 7-ray experiments. The black arrow 
indicates the systematic uncertainty on the energy scale. 



photon index of 0.15. Additionally, there is another sys- 
tematic uncertainty affecting the energy scale measurement 
at the level of 15-17%. Systematic uncertainties, related 
to the non-perfect matching of the observation conditions 
and the Monte Carlo simulations, have been estimated by 
studying the spectra reconstructed in different zenith an- 
gle ranges, analysis cuts, mis-pointing values, night sky 
background levels 1 10|. At all energies, the overall uncer- 
tainty is dominated by systematic errors rather than statis- 
tical ones. We fit the differential spectrum with a variable 
power-law: 

dN _^ / E X-^+^Logioli^) photons 
dEdtdA ^ ° V ITeV J TeVcm2s 

(1) 

where fo = (3.27 ± O.OS.tat) x 10^", a = 2.40 ± O.Ol.^t, 
and b = -0.15 ± O.Olstat- The residuals to the fit show sig- 
nificant deviations which are explained by systematic ef- 
fects ix^/ndf = 72/15). 

Figure [3] shows the spectral energy distribution (SED) for 
the same data points in Figure |2] and compares it to other 
measurements by lACTs. At energies above 10 TeV, given 
the large systematic error, it is not possible to disentangle 
the discrepancies between the results obtained by HESS 
ifTSl and those by HEGRA HI. 

In order to estimate the position of the IC peak we fit a 
variable power-law function to our data combined with the 
FermifLPsI data points taken from 1^. The fit accounts for 
the correlation in MAGIC data point and considers statisti- 
cal errors only. The obtained IC peak is 59 ± 6 GeV. 

4 The light curve 

In this section we present a time-resolved measurement, 
i.e. the light curve, of the 7-ray flux above 300 GeV from 
the Crab Nebula. This is meant to check the flux stabil- 
ity on timescales of days. Results are shown in Figure |4] 
The upper panel shows the light curve between October 1, 



2009 and March 31, 2011. In addition, the bottom panels 
are zoomed per-month portions of the bottom panel. The 
average flux of all nights F>3oogcV is: 

F>300GoV - (1.31±0.03,tat±0.173j,,)xl0-i"cm-2s-i. 

(2) 

The Crab Nebula flux is stable within the systematic un- 
certainty with a probability of 95%. The red shadow area 
in the panel marks the period of the enhanced 7-ray flux 
detected by both AGILE and Fenni in September 2010. 
MAGIC observed the Crab Nebula for one night during the 
same period, on September 20, 2010. The obtained flux is 
perfectly compatible with the average flux of all the nights 
within the statistical error. 

In addition to the sample used for the analysis presented up 
to here, we took some data between April 12 and 14, 201 1, 
for a total amount of about 140 minutes. These observa- 
tions were carried out under strong moonlight conditions 
with a special high voltage setting. A preliminary analysis 
of this second data set shows no increase in flux, as illus- 
trated in Figure |5] 
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Figure 5: Light curve of the Crab Nebula between April, 
12 and 14 during the flare detected at high-energies. 



5 Discussion 

The broad-band SED of the Crab Nebula has been matched 
with two different model calculations that are based on the 
model first suggested by [16J assuming a constant magnetic 
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Figure 4: Light curve of the Crab Nebula. 



field, B, as well as a solution of the ideal radial MHD flow 
||2]- The best match of the free model parameters (average 
B-field and magnetization of the flow at the shock, a) to the 
measured SED has been obtained using a rather weak field 
of S = (124 ± 6)AiG' lEl or a value a = 0.0045 ± 0.0003. 
Both model calculations are shown in Figure |6] together 
with the archival data as well as the new measurement pre- 
sented here. The spectral measurement obtained with the 
MAGIC telescopes covers for the first time the crucial en- 
ergy range at the peak of the IC component. None of the 
models provide a satisfactory match to the data if we con- 
sider statistical errors only. However, if taking systematic 
uncertainties into account, both models can fit the obser- 
vations. In case the systematic errors of the measurement 
can be further reduced, the new measurement can therefore 
help to lift a degeneracy and will require the introduction of 
additional structure in the B-field downstream of the shock. 
Specifically, the shape of the spectrum at « 100 GeV which 
is sensitive to the B-field in the region well outside the torus 
of the Crab Nebula. 

The constant B-field model offers the opportunity to cross 
calibrate the lACT measurements with Ferm;7LAT obser- 
vations. For the lACTs an energy scaling factor siact is 
introduced to correct the measured energy i?mcas to a com- 
mon energy scale E = sjact x -Emcas- This procedure 
reduces the uncertainty of the energy scale to that of the 
FermifLP^ 1 17|. The scaling factor is determined via a x^- 
minimization and is found to be siact = 1.050 ± 0.003. 
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